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Abstract 

Simple expressions are given for the Newtonian viscosity ?7Ar(</>) 
as well as the viscoelastic behavior of the viscosity rj((f),u)) of neutral 
monodisperse hard sphere colloidal suspensions as a function of volume 
fraction ef> and frequency u> over the entire fluid range, i.e., for volume 
fractions < <\> < 0.55. These expressions are based on an approxi- 
mate theory which considers the viscosity as composed as the sum of 
two relevant physical processes: rj((f>,(j) = f?oo(0) + Vcd(4',u), where 
"Hoo{4>) = VoX(4>) is the infinite frequency (or very short time) viscosity, 
with rj the solvent viscosity, x(4>) the equilibrium hard sphere radial 
distribution function at contact, and r) c d(4>,w) the contribution due 
to the diffusion of the colloidal particles out of cages formed by their 
neighbors, on the Peclet time scale rp, the dominant physical pro- 
cess in concentrated colloidal suspensions. The Newtonian viscosity 
Vn{4>) = *n{4>)U> = 0) agrees very well with the extensive experiments 
of Van der Werff et al and others. Also, the asymptotic behavior for 
large lo is of the form r)oa{4>) + A{4>){lotp)~ 1 / 2 , in agreement with these 
experiments, but the theoretical coefficient A{<f>) differs by a constant 
factor 2/x(4>) from the exact coefficient, computed from the Green- 
Kubo formula for r/((j),uj). This still enables us to predict for practical 
purposes the visco-elastic behavior of monodisperse spherical colloidal 
suspensions for all volume fractions by a simple time rescaling. 
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1 Introduction 



In a number of previous papers we have discussed the Newtonian viscosity as well as the 
visco-elastic behavior of concentrated colloidal suspensions, consisting of monodisperse neu- 
tral hard sphere particles' 1-4 ]. The motivation was to understand theoretically the very 
extensive viscosity measurements on colloidal suspensions carried out by Van der Werff et 
al in Utrecht' 5 ' 6 !. In particular, these experiments on carefully prepared systems seemed 
to be an ideal testing ground for the theory. In this paper a more complete and detailed 
account of the viscous behavior of colloidal suspensions over their fluid range will be given. 

Our theoretical approach is based on two physical processes related to the two widely 
separated basic time scales in a colloidal suspension: the Brownian time tb ~ 10 -8 s, dur- 
ing which a single Brownian particle forgets its initial velocity and the interaction time 
or Peclet time rp = a 2 /ADq ~ 10~ 3 s, during and beyond which Brownian particle inter- 
actions take place. Here a is the diameter of the hard sphere colloidal particles and Dq 
the Stokes-Einstein colloidal particle diffusion coefficient at infinite dilution. The viscosity 
is consequently considered as composed of a sum of contributions which take place on a 
short and a long time scale. Although the theory is constructed for concentrated colloidal 
suspensions with volume fractions 0.3 < <f> < 0.55, it appears that the theory also gives 
good numerical results for lower concentrations, so that effectively formulae are obtained 
which cover the entire fluid range < (f> < 0.55. Here <fi = where n is the number 

density of the hard sphere colloidal particles. 

The suspension is considered as a homogeneous fluid consisting of spherical particles 
immersed in a continuum solvent. As a consequence formulae derived for simple homoge- 
neous fluids in general - like the Irving-Kirkwood expression for the pressure tensor' 7 ' 8 ] or 
the Green-Kubo formula for the viscosity' 9 ] - are also assumed to be applicable here. The 
formulae for the viscous behavior are derived under a number of assumptions, which we will 
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try to justify physically as well as possible, but which, considering the complexity of this 
strongly interacting system, we have not been able to derive from first principles or justify 
completely. 

The two basic physical processes we referred to are 

1. at short times t < Tp <C tp and nonzero concentrations, the viscosity of the suspen- 
sion effectively increases when compared to that of the (pure) solvent viscosity 770 at infinite 
dilution, due to the finite probability to find two particles at contact; 

2. at long times f ~ rp > tb, the difficulty of a Brownian particle to diffuse out of 
the cage formed around it by its neighbors, characterized by a cage-diffusion coefficient 
D c (k;4>). 

As to 1., the probability to find two particles in the suspension at contact is given by 
the equilibrium radial distribution function at contact: g eq (a;4>) = x(<^)' 10 ', which follows 
from the canonical distribution of the hard sphere colloidal particles. As a result, the 
effective very high frequency viscosity of the suspension satisfies r/oo(0) = f/oX (</>)> a relation 
which is consistent with experiment over the entire fluid ranged . Similarly, the short time 
self-diffusion coefficient of the Brownian particles past each other is decreased from the 
Stokes-Einstein value Dq at infinite dilution, to a value D s {4>) = Do/x(((>), since x{<fi) a l so 
gives the increase in the binary collision frequency in a dense hard sphere gas in equilibrium 
as compared to that in a dilute gas. Also this relation has been confirmed by experiment W. 

As to 2., the cage diffusion coefficient D c (k;<fi) refers to the diffusion of a particle out 
of a cage formed by its neighbors when the particles are distributed periodically in the 
solvent with a wave number k. For concentrated suspensions one should bear in mind that 
a typical wave number is k « k* = 2ir/a, corresponding to a surface to surface distance of 
two neighboring Brownian particles of typically 1/10 of their diameter a, so that the particles 
"rattle" in their cages before they diffuse out in a time of the order of Tp ~ T c {k*]<p) = 
\/D c {k*;(f>)k* . In fig.l T c (k; <j>) / Tp is plotted as a function of k = ka for four values of (/). 
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T c (k;(f>) and Tp are clearly of the same order of magnitude, the pronounced maximum of 
T c (k; (f>) at k = k* corresponding to the "rattling in the cage". An explicit expression for the 
cage diffusion coefficient D c (k;(f>) has been obtained from kinetic theory! 11 !. Since D c (k;4>) 
also characterizes the decay of a spontaneous density fluctuation of wave number k in the 
suspension' 12 !, it can be measured by light or neutron scattering and the expression we give 
for it below has been shown to be in good agreement with such experiments' 13 ! . 

To incorporate the cage-diffusion process, i.e., D c (k;4>) into the theory, we need to go 
to a Fourier (i.e., k-) representation, while the starting point of our theory, the two particle 
Smoluchowski equation' 14 !, is expressed in ordinary (i.e., r-) space. This will introduce a 
fundamental difficulty in the development of the theory, since the impenetrability of two 
hard sphere particles, which is easily accounted for in r- space, will be violated in our theory 
in k-space, a point that will be discussed further below. 

The paper is constructed as follows. In section 2 we give the basic equations for the vis- 
cosity of the colloidal suspension and for the nonequilibrium pair distribution function of the 
colloidal particles to obtain this viscosity from a solution of the latter equation. In section 
3 this solution is used to obtain an explicit expression for the visco-elastic behavior r](<f), to) 
of the suspension. Section 4 gives a simple formula for the zero-frequency or Newtonian 
viscosity t]n(<P) = vi^i^ = 0)> while section 5 contains the visco-elastic behavior of the fluid 
for finite frequencies. In section 6 the approach of r]((f),uj) to its asymptotic value r/oo (</>), 
via a behavior ~ ^4(0)(o;rp) _1 / 2 , is discussed and exact results for the coefficient A((j)) are 
compared with our theory and with experiment. In section 7 the behavior of r](4>, to) for 
small oj is given and section 8 discusses a number of issues raised by the results obtained 
in the paper, especially in connection with the good agreement with experiment, in spite of 
the apparent neglect of hydrodynamic interactions between the Brownian particles. 
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2 Basic Equations 

The shear viscosity we are concerned with in this paper is defined as the linear response of 
the suspension to an applied shear rate j(t) = 7oe _lw * with finite frequency to and vanishing 
amplitude 70, or equivalently by 

Pxy{<t>, u, To, t) = -r](<l),u,j ,t)j(t) (1) 
Here P xy is the xy-component of the pressure tensor of the suspension, defined by 

P xy ((f),uj,jo,t) = P xy , s ((t>,7o,t) +P xy , d (4>,u, 7o, t) (2) 

where P xytS (4>, 70, t) is the static contribution {uj = 00) to the xy-component of the pressure 
tensor and P xy>( i{4>i UJ il0i't) the dynamic contribution given by^ 

1 N dV(ri) 

P xy ,d((p,UJ,"fo,t) = ~— < r O> Q r . tJ >n.e. (3) 

Here V is the volume of the system, rj the position of particle i (i = 1, N), = 
Yi — rj, V(rij) the interparticle potential between particles i and j at a distance = |ry| 
and the non-equilibrium average < > n . e . is taken with respect to a nonequilibrium distri- 
bution function derived from the iV-particle Smoluchowski equation for a suspension under 
shear rate Kinetic contributions to the pressure tensor are not considered in such a 

description of the system. 

The static contribution follows from the limit uj —> 00 when the dynamic contribution 
to the pressure tensor becomes zero, leaving in eq.(l) only 

P xy (<t),u = oo,7o,t) = P xy ,s{4>,lo,t) = -r?ooOh(*)- ( 4 ) 

Carrying out the implied integration on the right hand side (r.h.s.) of eq.(3) over the 
positions of all (N — 2) particles, but the particles 1 and 2, introducing center of mass and 
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relative coordinates by R = (ri + Y2)/2 and r = ri — r2, respectively, and carrying out the 
integration over R, one obtains for the dynamic contribution to the pressure tensor 

n 2 f dV(r) 
P xy ,d((f),uj,jo,t) = — — J drg(r;<f),uj,-f ,t)x — (5) 

This gives with eqs.(2) and (4) the following expression for the total pressure tensor: 

n 2 f dV(r) 
P xy ((f),uj,jQ,t) = -Voo(<f>h(t) ~ — J drg(r;(/),uJ,-fo,t)x — (6) 

Here n 2 g(r; (ft, uj, 70, t) is the nonequilibrium pair distribution function, giving the average 
number of colloidal particle pairs at a separation r in the suspension at a number density n 
of the colloidal particles, so that g(r;c/),uj,"f ,t) is the nonequilibrium generalization of the 
radial distribution function g eq {r;4>) in equilibrium, when 70 = 0. Introducing then: 

g(r;4>,u,j ,t) = g eq {r;4>) + Sg(r; </>, u, jo)e~ ilJt (7a) 

we have for 70 — > 

<^(r; 4>, u, 70) = 70<5ff(r; 4>, u) + O(7o) (7b) 

and one finds from eq.(6) that in the limit of vanishing shear rate 70 — ► 0, P xy {4>, w, 70, t) is 
proportional to since the contribution of g e q(r; <p) vanishes. Then in eq.(l), the viscosity 
rj((f>,(jj) = lim 7o ^o ? ?( < / ) ) ^1 70)0 is independent of 70 and t and given by: 

1 f dV(r) 
r]((f), u) = r/oo^) + -n 2 J dr5g(r; 4>, u)x — (8) 

An approximate equation for 5g(r; <f>, to) can be obtained in the following way. Neglect- 
ing the hydrodynamical interactions between the Brownian particles transmitted via the 
solvent, the iV-particle Smoluchowski equation for this case in a shear field j(t) can be inte- 
grated over the positions of all (TV — 2) particles but the two particles 1 and 2. This leads to 
an equation for the nonequilibrium pair distribution function, involving the nonequilibrium 
three particle distribution function. Neglecting the latter, i.e. restricting ourselves to low 



densities (to O(0 2 )), transforming to center of mass and relative coordinates of the two par- 
ticles 1 and 2, neglecting the dependence on the former, i.e., assuming spatial homogeneity 
and using g eq (r;(p) = exp(— f}V(r)), one obtains the following equation for g(r; (f>, u, 70, t): 

[- + 2/3 A, V • F(r) - 2D V 2 + l(t)x—} g(r; </>, iv, lo ,t) = (9) 

Here F(r) = — VV(r) is the force on particle 1 at a separation r from particle 2, (3 = 1/&bT, 
with ks Boltzmann's constant and T the absolute temperature. Eq.(9) has been considered 
for charged colloidal suspensions in the stationary state, i.e. for u = by Dhont et all 15 !. 
With eq.(7), eq.(9) can be written as an equation for 6g(r;(/>,(j): 

[ — iu! + 2pD V ■ F(r) - 2D V 2 ]5g(r; 0, w) = -x—e~ pv ^ (10) 

dy 

which has been solved exactly by Cichocki and Felderhofl 16 ! for hard sphere particles 
(cf. Appendix A). 

From now on we shall explicitly use a hard sphere potential unless specified otherwise. 
Neglecting then the force term on the left hand side (l.h.s.) of eq.(10) and taking the Fourier 
transform of eq.(10) with respect to r, an equation is obtained for: 

SS{k; <j>,w) =n J dre^dgir; 0, lo) (11a) 

Using that 

S eq (k; <f>) = l+nj dre ik - r [<? e(/ (r; 0) - 1] (lib) 

is the static structure factor in equilibrium in general, the equation for ^^(k; 0, u) derived 
from eq.(10) becomes: 

[-iu + 2D k 2 ]5S(k; 0, u) = 240^j 2 (ka) (12) 
where j2{ka) is the spherical Bessel function of order 2^ 17 \ 
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As pointed out in the Introduction, the neglect of the force term (which is only justified 
for r > a) in taking the Fourier transform of eq.(10), is the source of an error in the theory 
used in this paper to obtain the viscosity rj(4>, uo). A more detailed discussion of the nature 
of this error, its consequences and a way to partially correct for it can be found in section 
6 and Appendix A. 

Eq.(12) is only valid for dilute suspensions where g eq (r; 0) = exp(— (3V(r)), i.e. S eq (k; (j)) = 
1 — 24(f)ji(ka)/(ka), and where the basic diffusion process of the two particles is free diffu- 
sion, represented by the term 2Dok 2 on the l.h.s. of eq.(12). 

In order to obtain an equation for concentrated colloidal suspensions we make two 
corrections, a static one and a dynamic one. The first one replaces the low density expression 
for the equilibrium static structure factor, used to derive eq.(12) from eqs.(10) and (11), 
by the complete S eq (k;(f)) for concentrated colloidal suspensions. For the second correction 
we postulate that for such suspensions the basic diffusion process is cage diffusion rather 
than free diffusion. An expression for the relaxation time r c (fc; <f) for cage diffusion for 
concentrated colloidal suspensions has been derived before from the kinetic theory of a 
dense fluid of hard spheres, as the (scaled) reciprocal of the lowest eigenvalue D c (k; 4>)k 2 of 
a linear generalized kinetic operator, discussed elsewhere' 11 " 13 ! : 

D c {k,4>)k 2 = ,, d(fe) (13) 



T c (k;(f>) " ' ' x{4>)Seq{k\4>) 
Here D c (k; 4>) is the cage diffusion coefficient, S eq (k; <j>) is again the equilibrium static struc- 
ture for all 4> and d(k) = 1/(1 — jo(k) + 2j2{k)) a combination of spherical Bessel functions 
jt(k) of order t = and I = 2^. Writing 

d*) -"»**> (14) 

the frequency ujH{k;4>) is the half- width at half height of the dynamical structure factor 
S eq (k;uj) of the suspension in equilibrium, which is the quantity that can be measured in 
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light scattering experiments. The equality (14) is very well supported by experiment^ 13 !. 
Then eq.(12), becomes with eqs.(13) and (14): 

d 

[-no + 2u H (k; cj))]6S (k; (f>, to) = k y -^- S eq (k; 4>) (15) 
which has the solution: 

XQflr A ^ k xk y S^{k]4>) 

ss^M = — 2u)HM) _ iu (16) 

where S' eq (k; </>) = dS eq (k; 4>)/dk. 

We note that S eq (k; <f>) has a very sharp maximum at k ~ k* = 2ir/a at high densities' 13 ^ 
indicating a quasi periodic ordering of the colloidal particles on the length scale a in cages. 

Eq.(16) for SS(k; <f>,u) can be used to compute r}{(f>,u) with eqs.(8) and (11). This will 
be shown in the next section. 

3 General expression for the viscosity 

In order to use eq.(16) for ^(k; <j>,u) to compute r](4>,ui) we must Fourier transform eq.(8). 
For a hard sphere potential such a transformation is not possible. Therefore we replace 
in the spirit of the mean spherical approximation! 18 ], V(r) on the r.h.s. of eq.(8) by the 
equilibrium hard sphere direct correlation function C eq (r; 4>), i.e., 

V(r) - -k B TC eq (r; <j>) (17) 

As discussed in Section 6 and Appendix A, this replacement corrects partially for the neglect 
of the force term on the l.h.s. of eq.(10), which leads to unphysical contributions from 
overlapping particle configurations. Fourier transforming then eq.(8) by using Parcival's 
theorem on the r.h.s. and that the Fourier transform C eq (k;4>) of C eq (r;4>) is related to 
S eq (k;(f>) by: 

nC ^« =1 -^(b) <18) 
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one obtains straightforwardly from eqs.(8) and (11) the expression: 




6S(k;(t),uj) 



(19) 



Substituting eq.(16) into eq.(19) we obtain, after an angular integration in k-space: 



for the visco-elastic behavior of the suspension. 

In so far as the integrand in the second term on the r.h.s. of eq.(20) contains the 
eigenvalues (u;#(/c;0)) and amplitudes (S' eq (k;(f>)/S eq (k;(t))) of two cage diffusion modes, 
this term can be called a mode-mode coupling contribution to the viscosity. The difference 
with the usual mode-mode coupling contributions is that here two cage-diffusion modes, 
which describe the diffusion process in and out of two neighboring particles' cages, rather 
than two hydrodynamic modes (as occur in the long time tails or vortex diffusion^) are 
used. We also note that the same expression (20) for rj((j),u>) can be derived for ui = 0, 
by a direct application of mode-mode coupling theory to the Green-Kubo expression for 
7](<p,uj = 0)[ 20 K Since the complete derivation appears not to be in the literature, we briefly 
sketch it in Appendix B. For the concentrated suspensions we are mainly interested in here, 
the most important contributions to the integral in eq.(20) come from values of k « k*. 

We note that the /c-integral on the r.h.s. of eq.(20) is convergent for all w, since the 
integrand vanishes for k — > and the asymptotic behavior for k — > oo is ~ k~ 2 , as for large 




(20) 



s eq {k-4>) = i-i^ x {4>) 



ka 



[l + 0(fe- 2 )]; 



(21a) 



S' (fc;# = 24fo($ 



.72 

k 



[l + 0(k' 2 )}; 



(21b) 




k 2 [l + 0(k- 2 )] 



(21c) 
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This implies that the second term on the r.h.s. of eq.(20) vanishes for to — > oo, as it should, 
since r)((f), oo) = 77oo(0) by definition. 

We still have to obtain r?oo(^)) in order to compute r]((f),u). One often writes r?oo(0) = 
r/p 21 ^, i.e., equates r?oo(0) to the pure solvent viscosity, but this seems only correct for dilute 
solutions. For concentrated solutions, we propose to set: 

Voo(4>) = Vo X (<i>) (22) 

implying that the effective viscosity of the suspension at very high frequencies is not only 
determined by the solvent viscosity but increased by the fraction of colloidal particle pairs 
at contact, x(4>)- These touching, i.e. colliding particles, increase the effective viscosity 
proportional to the number of such pairs present in the suspension, because they increase 
the viscous dissipation in the suspension due to the instantaneous exchange of momentum 
during their collisions, no matter how short the time scale. They constitute therefore an 
instantaneous contribution to rj{(p,uj). Since! 10 ] 

X (<P) = 1 + \<t> + 4.590 2 + O(^) (23) 

eq.(22) reduces to the usual expression for rf^cj)) at small concentrations (see also section 
8, sub.3). 

In Fig.2 the behavior of r/ oo (0)/?7o = x(0) is compared with the reduced viscosity mea- 
surements by Van der Werff et al' 5 ' and Zhu et al' 22 ! at very high frequencies for <f> over the 
entire fluid range < </> < 0.55. Here we used the Carnahan-Starling approximation! 10 ] 



*<« ^ <24) 



which is very accurate for all such 4>. The agreement between theory and experiment is 
good, thus confirming eq.(22). We note, however, that a theoretical derivation of eq.(22) is 
lacking (see section 8, sub 3). 



11 



We also included in fig. 2 the values for t/oo(<^) as obtained by Cichocki and Felderhof ^ . 
These values differ slightly from those used by Van der Werff et al, since they obtained 
Voc(4>) by fitting the tails of the data for large u> to r?oo(</>) + A{4>)y/u)Tp, instead of using 
a fit for all uj. We used Cichocki and Felderhof 's values for r?oo(0) throughout the paper 
(cf. table 2). 

We remark that eq.(20) with eq.(22) and all the equations following from them, like 
eq.(25) in the next session, contain no adjustable parameters and are comletely determined 
by those charaterising the system: the viscosity of the solvent r/o, the volume fraction (j) (or 
equivalantly the number density n) and the diameter a of the colloidal particles. 

In the next two sections, we will compare the concentration dependence of the eq.(20) 
for the Newtonian viscosity t]n((P) = i]((f>,w = 0) and the concentration and frequency 
dependency of rj((f),oj) of eq.(20) with the experimental results of Van der Werff et al and 
others, respectively. 

4 Newtonian viscosity 

Setting oo = in eq.(20) and using eqs.(13), (14) and (22), we obtain the following simple 
expression for the Newtonian viscosity: 

K,(^)] 2 

where k = ka and the Stokes-Einstein relation 



toW=toW[1+ j_/;^k^i ] (25) 



D„ = (26) 

has been used. 

Although the expression (25) for t]n{4>) has been derived for large <p (0.3 < <j> < 0.55), 
where cage diffusion is the dominant finite time contribution to the viscosity (via eqs.(13) 
and (14)), eq.(25) nevertheless appears to describe the ^-dependence of r/N((j)) for small and 
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intermediate concentrations also, due to the presence of the ?7ox(0) term (cf.fig.3). Fig.3 
also shows that the cage diffusion describes the very rapid increase of r/jv(0) with </> for 
0.40 < (f> < 0.55 very well. 

Eq.(25) has been evaluated using the Henderson-Grundke correction! 24 ] to the Percus- 
Yevick equation for the computation of the hard sphere S eq (k; </>) and S' eq (k; (j)). A convenient 
Pade approximation of 7?at(0) for practical use for all < (f> < 0.55 is: 

within a relative accuracy of less than 0.25%. This approximation yields for r]M{4>) t ne 
correct Einstein coefficient |^ as well as the same coefficient of 0((f> 2 ) as eq.(25). 

Cichocki and Felderhof have obtained on the basis of the pair Smoluchowski equation 
exact results for rj(4>,u>) to 0{4> 2 ). Their result to 0{4> 2 ) for r]N((f>) is, without Brownian 
motion contributions' 25 ! : 

VN (4>) = l + ^ + 5.OO0 2 (28a) 
while, with Brownian motion contributions they find[ 26 h 

VN {4>) = l + ^ + 5.910 2 (28b) 

This can be compared with the approximate result we obtain from eq.(19): 

WW = 1 + |^ + 6.O30 2 (28c) 

where the term 6.03^ 2 contains a contribution 4.59</> 2 from ^oo(^) and a contribution 1.440 2 
from the second (mode-mode coupling) term between the square brackets on the r.h.s. of 
eq.(25). Since for (f> < 0.25 the cage-diffusion contribution to rf(4>;u) can be neglected, 
eq.(22) then reduces to r]N(<fi) = r/oo(0) = i]ox(<fi)- The eqs.(28b) and (28c) both give then 
a very good representation of the experimental values for 7/Ar(</>). 
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5 Visco-elastic behavior 



For r]{4>,uj) of eq.(20) is complex, so that the visco-elastic behavior of the suspension 

can be written in the form: 

ti&u) = r/(<f>,u;)+ir/'((f>,<j) (29) 

where rj'((p,uj), rj"{<p,uj) are the real and imaginary parts of rj((f),uj), respectively. It is 
convenient and customary' 5 ] to consider, instead of rj'((p,uj) and rj"((f>,uj) reduced quantities 
defined by: 

* / , \ rfi^u) -i]((p,oo) _ r/(<f>,u) - ?? oo (0) 

Vr{ ^ w) " ^,0-^,00) " V n(4>)-Voo(<P) (30a) 

and 

W,^) = — 7-t 7-7 (30b 

where the reduced real part r]* R ((j),uj) varies as a function of uj between 1 (for to — > 0) and 
(for u> — > 00) for all </> and rjj(4>,uj) vanishes for uj — > and — > 00, exhibiting a maximum 
in between. In fig.4, rf R ((f), uj) and ry|(0, u;) are compared with the experimental data of Van 
der Werff et al, as a function of a reduced u> for all available 4> for 0.44 < 4> < 0.57' 5 ' . As Van 
der Werff et al state, the values they find for the reduced quantities rj R (<f>,u) and r)j(cf),uj) 
are very weakly dependent on (f>, which is consistent with the crowding of all experimental 
points around the theoretical curves, inside the experimental errors. The scaling of uj for 
the experimental data was performed in the same way as was done by Van der Werff et al 
by fitting the data for large uj to the expression (cf. section 6): 

Vit(<l>,u)=vK<l>> u ) = ^r^ , 1 7 ,x ( 31 ) 

2TT ^UJTl{(t>) 

where t\{4>) is a phenomenological time for the experiments. The T\{4>) used for the theo- 
retical results is given in section 6, eq.(33). 

Nevertheless a more detailed comparison of rj RI ((j),u>) as a function of <\> can be made, 
although the large experimental uncertainties of the data and the difference in the basic 
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inputs in the theory ((/> and 770) and experiment (a, c and r/o, with c the weight concentra- 
tion of the colloidal particles) complicate considerably a compelling detailed comparison of 
theory and experiment. Examples are given in fig. 5. In the same figure the results of a 
general phenomenological description of the visco-elastic behavior of colloidal suspensions 
due to Cichocki and Felderhof are given I 23 ] . This description is based on a three pole ap- 
proximation in the complex y^-plane, whose location is derived from the experimentally 
measured values i] e ^ p ((/)), r]^ p (4>) and three additional parameters, one of them being a re- 
laxation time. From these three poles the ?/(</>, to) and r/'((/>, w) as a function of uj can be 
derived. For the three concentrations <ft = 0.44, 0.46 and 0.53, for which their procedure 
could be implemented, 7]'((j),uj) and rj"((p,uj) are consistent with our results within the ex- 
perimental errors. As was shown by Cichocki and Felderhof, the strongly deviating cloud 
of points near uit\{4>) w 1 in the imaginary part of the reduced viscosity ri}(4>,u;) (cf.fig.4b) 
can be disgarded, since they violate the Kramers-Kronig relations between the real and the 
imaginary part of rj{(p,uj) and must therefore be erroneous! 23 ]. 

6 Large u;-behavior 

For large u, eq.(20) for rj(4>,uj) can be written as: 

r?(0, u) = Vqo (<P) + % 2 X 5/ % -7^= (1 + i) + O(-) (32) 

where the square root singularity for uj — > 00 is induced by the large fe-behavior of the 
integrand on the r.h.s. of eq.(20), as given by eq.(21). We note that the correction 0(— ) 
is an exact result for low concentrations to 0{4> 2 ) (cf. Appendix A) and is consistent with 
what is found in the mode-mode coupling approximation. 

Using eq.(32) in eq.(30) and comparing with eq.(31) gives for T\{<j>) the theoretical 
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expression: 

which is plotted in fig.6 and is consistent with the experimentally used n(0) up to about 
4> ~ 0.55, averaging at a value of about rp/4 (cf. section IV. B in ref.5). The systematically 
too low theoretical value of t\ ((f)) corresponds to the systematically too high theoretical 
value of the coefficient of the w _1//2 -singularity in eq.(32) as compared with the exact value 
given in eq.(41) below. 

In fact, in order to investigate this behavior further, an independent evaluation of rj((j), u) 
for large oj was made, starting from a Green-Kubo like formula for r]((j),uj) rather than from 
eq.(8): 

V (<I>,lj) = V oo(<t>) + / dtp v (t- fie"* (34) 
v Jo 

Here the stress-stress auto correlation function p v (t) is defined by: 

p v (t;cf>)=<^ xy e nt ^ xy > eq (35) 

where the brackets denote an equilibrium ensemble average. Here, instead of using the 
microscopic pressure tensor (the expression within the square brackets of eq.(3) in section 
2), we use the in this context more customary microscopic stress tensor E?^, which is equal 
but opposite in sign and can be written as: 

N 

>■!■„ >:..,. !■), , (36) 
i=i 

with Fi = -Vi®(r N ) the force on particle i (V; = <9/<9r;), $(r N ) = Ei<j=i v ( r ij) the total 
potential energy of the colloidal particles and 

N 

n = D s J2[Vi + PFi\-Vi (37) 
i=i 

the A^-particle Smoluchowski operator! 27 ] with Do replaced by the short time self-diffusion 
coefficient D s (4>) to make eq.(34) applicable to all fluid densities. This is further discussed 
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below. For N = 2 and x{4>) = 1 the adjoint operator occurs in the pair Smoluchowski 
equation (eq.(9)). 

The short time behavior of p v (t; <p) determines the large oj behavior of rj(4>, u). Since for 
hard spheres the interparticle potential is singular, one determines the short time behavior 
of p v (t; <p) by first using a soft potential Vi{r) = e(^) e , where e is the two particle interaction 
energy for r = a, and then letting I — > oo, so that V((r) approaches a potential between two 
hard spheres of diameter a. For £ — ► oo, one can then derive for p v (t, 4>) the expression! 28 ] : 



. . 2Tm 2 Va 3 x(<P)t , 
Pv(t; <t>) = 15/j2 r(t*) (38a) 



with 



•(**) = I" '<fce-VK' a & + *-' 2 >& a (38b) 
Jo 



where 

2D tf 2 

t* = (38c) 

The leading term of r(t*) for \im t ^ lim^ooi i.e., t* ~ tl 2 — ► oo, which determines the short 
time behavior of p v (t;cj)) for a hard-sphere potential, was obtained by M. J. Feigenbaum 
and reads' 28 ': 

KO = (39) 

V 7Tt 

Using eqs.(34), (38) and (39) and the Stokes-Einstein relation (26), one obtains for rj((/),u) 
for large u and for a hard sphere potential for all <j> the exact expression: 

D, 



18 

rj((/),u) ~ ^(^ + y0 2 x(<A)% 



*/ 2 1 + i 

^? <40) 



Using then that D s {4>) = Dq/x{4>) (cf. section 8, sub. 3) one has: 

18 1 

W)~ 7/^(0) + _0 2 X (^3/2 r?o __ ( l +i) (41) 
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Eqs.(32) and (41) are both compared with the experimental data for large ui and for most 
experimental values of 4> in fig. 7. We emphasize that in order to get agreement with exper- 
iment it is necessary to replace the low density Stokes-Einstein diffusion coefficient Do by 
the self-diffusion coefficent D s (<fi) in the basic Smoluchowski operator (cf.eq.(37) and fig. 7). 
We also emphasize that the exact result of eq.(41) constitutes a generalization of Cichocki 
and Felderhof 's low concentration result to all concentrations in the fluid range. A detailed 
derivation of eq.(41) will be given elsewhere' 28 ' . 

It is clear that the experiments agree very well with eq.(41) and not with eq.(32), 
consistent with the systematically lower theoretical values of ti{4>) in fig. 6. This could 
well be related to the approximations made to obtain eq.(32): (1) the use of the complete 
S eq {k;4>) (i.e. for all (f>) in the two particle eq.(15) and the use of ujH{k;4>) as the only 
basic relaxation time; (2) the replacement of the potential V(r) in eq.(8) by the direct 
correlation function C eq (r;(f)) and (3) the neglect of the force term on the l.h.s. of eq.(10) 
and consequenly the correct boundary condition of hard sphere impenetrability incurred by 
the Fourier transform from eq.(10) to eq.(12) (cf. Appendix A). 

The first approximation was intended to incorporate in the calculation of r]((/),uj) contri- 
butions due to more than two isolated particles, i.e., correcting for the neglect of the three 
particle distribution function in the eq.(9) for g(r; (f), uj, 70, t). 

As pointed out before, the second approximation, is necessary to perform a Fourier 
transform of eq.(8). It also corrects partly for the unphysical contributions from overlapping 
particle configurations, due to the neglect of the proper hard sphere boundary condition 
(cf. Appendix A). We remark that the Fourier transform of eq.(8) was due to the necessity 
of introducing the relaxation times r c (/c; </>) related to the cage diffusion for concentrated 
colloidal suspensions, which have only been determined for periodic particle arrangements, 
characterized by a wave number k. However, neither of these two approximations seem to 
be responsible for the incorrect asymptotic oj-behavior of rj{4>,uj). 
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As to the third approximation, if we compare eqs.(32) for low densities, i.e. x(0) = 1> 
with the exact solution for rj(4>,Lv) obtained by Cichocki and Felderhoft 16 ] to 0(4> 2 ), we see 
that the second term on the r.h.s. of eq.(32) is smaller by a factor 2. Cichocki and Felderhof 
considered eq.(10) with the correct hard sphere boundary condition in r-space and solved 
it exactly for all t. If we solve eq.(10) in the same manner but neglect the force term on the 
l.h.s. (cf.Appendix A), we obtain, however, eq.(32) in the limit of large u with = 1. 
This suggests that the third approximation, the neglect of the force term on the l.h.s. of 
eq.(10) and the ensuing violation of the proper hard sphere boundary condition in real space 
in making the Fourier transform from eq.(10) to eq.(12) is the main reason for the erroneous 
expression (32). 

We note that the eqs.(32) and (41) show that the difference between the exact and the 
mode coupling result for the coefficient of u^ 1 / 2 is a constant factor 2/x{4>)- This only 
affects the approach to uj = oo, not ??oo(</ ) ) itself, and is of no influence if one plots the mode 
coupling theory on the phenomenological time-scale loti{4>) using eq.(33) (cf.fig.5). This 
may be of practical importance for predicting the visco-elastic behavior of concentrated 
colloidal suspensions since the scaling in time does not affect the Newtonian behavior of the 
viscosity t 29 l. 

7 Small a>behavior 

For low densities to 0((f> 2 ) the small u, or long time, behavior of r/(4>,uj) follows from 
eqs.(20), (21) and (29) to be: 

7/o 25 175 

r)"{d),oj) 48 . . , 9 

■L^-L = —(utp)4> 2 + .... (42b) 
7/o 175 
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This can be compared with the exact results of Cichocki and Felderhof' 16 ' to 0(4> 2 ) for 
uj -> 0: 

= { ^ _ ™ {u)Tp?}(t ? + .... (43a) 
V"{<t>,u) 8 2 

— = — </» (wt p ) + .... (43b) 

'/O 10 

The agreement of eqs.(42a,b) with eqs.(43a,b) for small uj and low concentrations, in partic- 
ular of the coefficient of {ujtp) 2 in the real parts, is better than that of eq.(32) and eq.(41) 
for large uj. This is probably due to the fact that the neglect of the proper hard sphere 
boundary condition in the mode-mode coupling theory is more serious for a description of 
the short time behavior than the long time behavior of the suspension. We remark however 
that the difference in the first terms on the r.h.s. of the eqs.(42a) and (43a), i.e. 36/25 
and 12/5, respectively, is a direct consequence of the violation of the proper hard sphere 
boundary condition (cf.Appendix A, in particular eq.(A.25)) 

8 Discussion 

1. The w-dependence of rj((f),uj) is well represented by eq.(20) for all (f> on the phenomeno- 
logical time-scale T\(4>) or if plotted as a function of ujtp, when an over-all shift to the 
theoretical curves of 2/x{4>) is applied! 29 ]. This is due to the fact that the asymptotic 
mode-mode coupling result (32) for the large uj behavior of r)((f), uj) is not correct, because 
of the incomplete incorporation of the hard sphere impenetrability in the theory. The mode- 
mode coupling contribution to rj((f),uj) should be best for values of uj around ujT\((p) « 1, 
where there are rather few experimental points. It would be interesting therefore if a more 
detailed comparison between theory and experiment could be made in this u;-regime, to 
obtain a more appropriate test for the validity of the mode-mode coupling theory used 
here. 
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2. The result (20) for rj((f>, to) is based exclusively on the instantaneous time behavior 
of 77oo(</>) and the cage-diffusion relaxation mechanism. From the agreement of rj(4>,u)) and 
Vn{4>) with experiment, it would seem that these two physical processes essentially suffice to 
understand the Newtonian as well as the viscous behavior in the entire fluid range of hard 
sphere colloidal suspensions. That this agreement occurs without considering explicitly any 
hydrodynamical interactions between the colloidal particles in the theory presented here 
may appear rather puzzling. We do not have an explanation for this, other than that at 
high concentrations, where 0.3 < (f> < 0.55, the surface to surface distance between the hard 
spheres is so small, that a "quenching" of hydrodynamical effects is not unthinkable. 

3. There may, however, be a deeper justification for the neglect of the usual hydrody- 
namical interactions in our theory. It seems that in a number of cases the same dependence 
of a physical quantity of the suspension can be obtained by theories with and without hy- 
drodynamical interactions between the Brownian particles. In this respect the following 
two observations are relevant. 

(a) The concentration dependence of the infinite frequency viscosity r?oo(0) as well as of 
the Newtonian viscosity 7?Ar(</>) for low and intermediate concentrations < (f> < 0.45 are 
described by our relations (cf.eqs.(22) and (25)): 



?7oo(0) 



mx{4>) = 

r ?o [l + ^ + 4.590 2 + O(0 3 )] 



(44a) 



and 




= % [l + ^ + 6.O3^ 2 + O(0 3 )] 



(44b) 



respectively. The r.h.s. of eqs.(44a) and (44b) can be compared with Beenakker's expression 



[30]. 
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rf"(0) = lim k ^ [r l (k; ( j>)} = 

= 77 O [l + ^ + 4.840 2 + O(<A 3 )] (44c) 

for, what he calls, the effective viscosity. Beenakker's rf$ f (4>) is derived from a wave vector 
dependent viscosity rj(k;(f)), a complicated functions of k, by using the quasi-static Stokes 
equation to describe the motion of the fluid, neglecting inertial effects. This implies, as he 
pointed out, that his equation is valid for Tp < t < Tp. Our relations (44a) and (44b), 
however, are valid for t < Tp and t > Tp, respectively. Thus his result (eq.(44c)) can be 
regarded as between eq.(44a) and eq.(44b) (cf.fig 8a). While for low concentrations the 
difference between the three expressions (as well as eqs.(28a) and (28b)) is marginal, since 
it does not appear to be relevant for comparison with experiment, we emphasize that the 
strong experimental increase of the Newtonian viscosity for higher concentrations <j) > 0-3, 
can only be described by the integral on the r.h.s. of eq.(44b) (cf.figs.2 and 8a). 

(b) Also, the concentration dependence of the short time self-diffusion coefficient D s (<j)) 
for low and intermediate concentrations < < 0.45 can be equally well described, within 
the experimental uncertainties, by our relation: 

D. W - J* (45a) 

as by the Beenakker and Mazur expression! 31 ! : 

D s {4>) =]im k ^ 00 D(k;<i>) (45b) 

where D(k; <f>) is a wave vector dependent collective diffusion coefficient, which is, like 
r](k; <f>), a complicated function of k. While our relation (45a) for D s ((f>) is valid for t < Tp, 
Beenakker and Mazur's expression (45b) is, like for their viscosity, valid for rp < t < rp. 
On this larger time-scale D s (4>) will contain extra, in his case, hydrodynamic contributions 
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in addition to our instantaneous contributions, leading to slightly larger values for the short 
time self-diffusion coefficient. The same obtains for the experiments of van Megen et al' 45 ! 
and Pusey and van Megen t 46 l (cf.fig 8b). 

Beenakker and Mazur consider only purely hydrodynamic interactions between the par- 
ticles, in that they study the hydrodynamical effect of a number of stationary particles on 
the motion of one moving particle. In our case no hydrodynamics enters explicitly at all, 
essentially only molecular considerations are used. For short times the (static) equilibrium 
radial distribution at contact %(</>), derived from the canonical distribution of the colloidal 
particles in equilibrium, occurs, yet a comparable agreement with experiment is obtained. 
For long times there is an extra (dynamic) contribution due to the increasing difficulty for 
a particle to diffuse out of the cage formed by its neighbors. 

(c) We believe that for a complex system like a colloidal suspension there could be ap- 
parently very different alternate descriptions of the same phenomena. Perhaps the simplest 
and most striking example of this is the observation that Einstein's low concentration result 
for the viscosity of a colloidal suspension, derived from Stokes hydrodynamics ' 32 1 

^ = l + h + Otf) (46a) 
Vo * 

can also be obtained, using an Einstein relation (cf.eqs.(44a) and (45a)): 

— =dTS = 1 + ^ + ^» < 46b > 

770 D s ((p) 2 

Although these equivalent alternate descriptions of colloidal suspension properties - and 
especially eq.(46b) - could well be a fluke, a deeper origin cannot be ruled out in our 
opinion either. 

In fact, for the equivalence of Einstein's expression (46a) and our (46b) the following 
physical argument can be given. 

Felderhof has shown' 33 ] - and it also follows from the Green-Kubo expression (34) - that 
rj((j), uj) = 770 [1 + §</> + T]2{uj)<fi 2 ]. Therefore the first two terms in the expansion of r](4>, u) in 
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powers of <p are independent of uj. This implies that when computed for any uj they should 
give the same answer: r/o[l + §</>]■ 

Einstein - as represented in Landau-Lifshitzt 34 ! - did the computation for uj = 0, i.e., he 
used a long time stationary state hydrodynamic calculation, to obtain the extra resistance 
of the suspension to shear from the change of the velocity field of the fluid due to a single 
Stokesian hard sphere particle placed in it. 

We propose to do a computation at uj = oo, i.e., for a very short (in fact, instantaneous) 
time. Then the placing of one particle - or even many mutually separated particles - in 
the solvent will not have any effect on the viscous resistance of the suspension. The only 
way the presence of the particles can produce an extra flow resistance is from pairs of 
particles (already) in contact, where an "instantaneously" collision takes place adding to 
the viscous dissipation in the suspension. Therefore, for uj = oo the increase in the effective 
fluid viscosity as a function of 4> will be given by the relative increase in the number of 
particle pairs at contact in equilibrium as a function of 4>, which is x(4>)- On the basis of 
this argument one would conjecture that for uj = oo, the increase in suspension viscosity, 
when compared with that of the pure solvent, would be x(^) f° r an 0) n °t just 1 + |(/> 
to 0((p). This conjecture is consistent with experiment (as shown in fig. 2) and should be 
derivable from kinetic theory! 35 ] . 

4. We also remark that the Einstein relation 

Do = ^~ (47a) 

appears to hold not only for infinitely dilute suspensions, but for all concentrations in the 
formW : 

XW) 37T7? oo (0)cr 

as can be seen in fig.8c. The physical reason for this seems to be that as long as the times of 
observation are sufficiently short (or the frequencies sufficiently high) , so that no significant 
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motion of the colloidal particles can take place, no hydrodynamical effects will occur, and 
only the instantaneous effect due to particles at contact - which does not require any time 
to occur -, i.e., x(0) wm be relevant. 

5. Recently Bradyt 36 ! has published a different model for the Newtonian as well as the 
frequency dependent viscosity. His results can be obtained from the low density result of 
Cichocki and Felderhof' 1 ^ (cf.Appendix A) with only two modifications: (1) a scaling of 
their exact solution (eqs.(A.2) and (A. 6)) for the low density two particle Smoluchowski eq. 
(10) (eq.(A.l)), by replacing the Stokes-Einstein diffusion coefficient Dq by the short time 
selfdiffusion coefficient D s ((fi) and (2) the addition of a factor g e q(r = a;4>) = x{4>) to the 
low density expression for the potential contribution of the viscosity in terms of the pair 
distribution function (cf. the second term on the r.h.s. of eq.(8)). This leads directly to 
Brady's expression for r)((f),uj) (cf.eq.(A.ll), which in our notation reads: 

r)((f>,u) = TfcoO) + T]o4> 2 a v {u))g eq {a; <f>) (48) 

Ds{<P) 

which reduces for uj = (with eq.(A.12)) to his expression for the Newtonian viscosity 

w(# 

12 D 

w(<t>) = »7oo(</>) + —r]o<t> 2 geq(<r; 4>) n ° (49) 
5 D s ((p) 

However, in his calculations Brady determines the three basic ingredients of his theory em- 
pirically: r/oo(0) is derived from measurements and Stokesian dynamics! 37 ' 38 ], while g eq (a; 0) 
is taken to be given by the Carnahan-Starling approximation eq.(24) for < 4> < 0.5 and by 
1.2(1 — ^/</> m ) -1 for (ft > 0.5, as derived from dense h.s. fluid computer simulations, where 
<j) m = 0.63 is the volume fraction of random close packing of hard spheres. Furthermore the 
relative short time self diffusion coefficient D s {4>)/Dq is taken from Ladd's computer simula- 
tions for < 4> < 0.45I 38 ! and from Phung's Stokesian dynamics simulations for cf> > 0.45t 39 ]. 
This leads to a curve for t]n (</>), as given by eq.(49), which is virtually indistinguishable 
from our t]n{4>) based on eq.(25) for < 4> < 0.55. We remark that eq.(49), with the just 
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mentioned determination of 7700 (</>), g eq (a;(p) and D s (4>)/Do, also describes very well the 
experimental data for tjn (<^)[ 5 ' 6 > 40 > 41 ] for 0.55 < <f> < 0.60, where the precise thermodynamic 
state of the suspension is not clear, while eq.(25) gives then too low values for r/^ (</>). Vir- 
tually the same result as Brady's description of 7?jv(</>) for < <p < 0.60 can be obtained by 
using in his eq.(49) for all <j>, our eqs.(22) and (45a) for r?oo(^) and D s (4>)/Do, respectively, 
as well as his representation of g eq (a;4>). It is clear that the precipitous increase of t]m{4>) 
for <fi > 0.55 is then a direct consequence of the pole in g eq (&; 4>) at <p = (f> m . 

However, for the visco-elastic behavior, when plotted as a function of ojt\ ((f)), Brady's 
results do not agree as well with the experiments of Van der Werff et al^ 36,42 ' . This may well 
be related to the fact that the basic ingredient of Brady's theory that causes the increase 
of 77iv(0) for large 4> is a static one, related to the behavior of g eq (&; 4>) ~ (1 — 4>/4>m)~ 1 as 
random close packing is approached, while in our theory it is a dynamic one: the increasing 
difficulty of diffusion of a particle out of the cage formed by its neighbors. It appears that 
only the latter one is able to account for the frequency behavior of r/(0; uj). The underlying 
physics of the two processes is therefore very different: while we use the typical high density 
mechanism of cage diffusion, Brady upgrades the low density physics by effectively scaling 
with g eq (a;4>) and D s (4>). 

We note that essentially the same mode- mode coupling term as in eq.(25) gives the steep 
viscosity rise at high densities for atomic liquids, since the atoms - like the colloidal particles 
- find themselves in cages, out of which they can only escape with increasing difficulty with 
increasing density^ 19,20 ] . 
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Appendix A 

Here we compare for low densities <j> — ► and hard spheres the exact dynamic viscosity 
r]((f>,uj) as obtained from eqs.(8) and (10) by Cichocki and Felderhofl 16 ! with the mode-mode 
coupling approximation rj mc ((j),uj) given by eq.(20). 

We first give the exact solution of eq.(10) for 5g(r;cf),uj) as obtained by Cichocki and 
Felderhof. For 0^0, g eq (r;(j)) = exp(— (3V(r)), so that eq.(10) reads: 

[-iw + 2A)V • {/3F(r) - V}]Sg(r; 0, u) = 0°^V {r)e^ v ^ (A.l) 

with V'(r) = dV(r)/dr. The solution of eq.(A.l) can be written as 

6g(r;<l>,u) = ^f(-;u)e-f lv W (A.2) 

Substitution of (A.2) into (A.l) and using that 

{/3F(r) - Vje-^M = (A.3) 

one obtains in the hard sphere limit V(r) = lim^^ V\ (r) = lim^oo e(r/a) 1 the following 
equation for f(u;u>), with u = r/a, 



with the boundary condition 



/'(i;-) (A.5) 



where f'(u;uj) = df(u;uj)/du. This boundary condition ensures that the r.h.s. of (A.l), 
which diverges sX r = a for hard spheres, cancels exactly a similar divergent term arising 
from F(r) on the l.h.s. The solution of (A. 4) with (A.5) is, for r > a(u > 1). 

,, x cr 2 fc 2 (au) 
219o afc 2 (a) 
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with k2(x) the modified spherical Bessel function' 17 ] of the third kind, 

k 2 ( x ) = e' x { X - 1 + 3x~ 2 + 3aT 3 } (A.7) 

and 

a = a(w) = (1 - (A.8) 

We note that for hard spheres f(r/a;uj) is continuous at r = a so that 6g(r;<j>,u) in (A. 2) 
shows a jump at r = a due to the factor exp(— (3V(r)) = 0(r — a) with 0(x) the Heaviside 
step function. In particular, 5g(r; (p,uj) = for r < a, reflecting the inpenetrability of two 
hard spheres. Next we substitute (A. 2) for 5g(r; (p,uj) in eq.(8) for rj((j);uj). Using that for 
hard spheres 

V \r)e~ pvir) = -k B TS(r - a) (A.9) 

one obtains straightforwardly 

2tt 

V (<j>; uj) = Voo(4>) ~ —k B Tn 2 a 3 f(l; u) (A.10) 
15 

Substitution of (A. 6) and (A.7) leads to the final result for (f> — > 0, 

f]((f); lo) = rjooic/)) + r]Q(f) 2 a v {uj) (A.ll) 

with 

36 a 2 + 3a + 3 
ay w = - — ^ —5 — (A.12) 

v ' 5 a 3 + 4a 2 + 9a + 9 v ' 

and a = a(w) given by (A.8). 

In the mode-mode coupling theory on the other hand, one neglects the force F(r) on 
the l.h.s. of (A.l), so that 5g mc (r; 0, u) satifies 

l-iiu - 2D V 2 ]5g mc (r; M = (3^V {r)e~P v ^ (A.13) 

The solution of this equation can be written in the form: 

xy v 

5g mc (r;<j),uj) = ^/ mc (-;w) (A.14) 
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Substitution of (A. 14) into (A. 13) yields the following equation for f mc {u;u): 
with boundary condition (e — ► 0) 

2 

f mc {l + e-u) - f mc {\ -e;u) = ^ 5 - (A.16) 

which follows from the r.h.s. of (A. 13) in the hard sphere limit, using (A. 9). Thus, 
f mc (r/a; u) is continuous for all r with a jump in its derivative at r = a given by (A.16). 
The solution of (A. 15) and (A.16) is for u < 1, 

f t \ ° 2 -K{a) i 2 {au) 

fmcWUj) = —— A.17) 
2Dq 1 + K (a) m 2 {a) 

and for u > 1, 

a 2 1 ko(au) 
fmc(u;u;) = — - -^-L A.18 

2Dq 1 + K (a) ak 2 {a) 

where a = a(u>) is defined in (A. 8), k 2 (x) in (A. 7), i2(x) is the modified spherical Bessel 
function of the second kind^ 17 ! , 

3 1 3 
hix) = (— , H — )sinhx ^-coshx (A. 19) 

x J x x 2 



and 

*(„) = -MfM (A.20) 

Thus (5<7 mc (r; 0, w) given by (A. 14), (A.17) and (A.18) is continuous for all r and nonvanish- 
ing for r < a, allowing two spheres to overlap. To exclude such unphysical configurations 
in eq.(8) for the viscosity r]((f),u) we replace V(r) by —kBTC eq (r; </>) (eq.(17)). Using that 
C eg (r;(f)) = exp(— (5V(r)) — 1 for <p — ► 0, dV(r)/dx in eq.(8) is then replaced by 

dV(r) _ ^ mr) dV(r) 



e-™^- (A.21) 



dx dx 
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and the factor exp(— (3V{r)) so obtained excludes the unphysical contributions in 5g mc (r; (ft, uj) 
for r < and thus partially compensates for the error made in the boundary condition of 
eq.(A.13) as far as -q((ft,uj) is concerned. Substitution of (A. 14) in eq.(8) with the replace- 
ment (A. 21) and using (A. 9) leads to 

2tt 

Vmc((ft,u) = Voc((ft) ~ —k B Tn 2 a 3 f mc {l;uj) (A.22) 
15 

completely similar to (A. 10) for r]((ft,uj). Using (A. 18) for f mc (l;uj) yields the final result 

rj mc ((ft,uj) = r/ooO) +T]Q(ft 2 av{u)) — (A.23) 

1 + K [a) 

with ay(w) given by (A. 12), K(a) by (A.20) and a = a(u) by (A.8). 

The result (A.23) for rj mc ((ft,u>) follows from eq.(20) provided one uses there the low 
density expression for S eq (k; (ft) and ujH(k) = D$k 2 . To compare the exact expression (A.ll) 
for r/((ft>,uj) with (A.23) for r] mc ((ft,uj) we note that for large frequencies lo — > oo,a — > oo 
(cf.(A.8)) and K(oo) = 1 (cf.(A.7), (A.19) and (A.20)), so that then 

rj mc ((ft,uj) - r]oo(<ft) = ^(rj{(ft,uj) - r?oo(0)) (A.24) 

For u -> 0, a -> (cf. (A.8)) and K(0) = 2/3, so that then 

3 

r] mc ((ft,uj) - r/oo^) = -(77(^,0;) - f]oo((ft)) (A.25) 

Thus it appears that the mode coupling theory underestimates the two particle Smolu- 
chowski contribution to rj((ft,Lu) by a factor 2 at high frequencies and 5/3 at low frequencies. 
The relevance of these factors is limited in practice since for low concentrations the main 
contribution to r)((ft,u) comes from ?7oo (</>)• For high concentrations the factor 2 is reduced 
by a factor x(0), due to the replacement of Dq by D s {(ft) in the two particle Smoluchowski 
equation (6) (cf.section 6). 
Appendix B 
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Here we derive eq.(20) for rj((f), u) directly, using the mode-mode coupling approximation 
(mmca) for concentrated suspensions 0.3 < 4> < 0.55, in analogy with what is done for 



of a given dynamical variable decay predominantly into pairs of modes associated with 
conserved single-particle or collective dynamical variables I 43 ] . If we restrict ourselves to the 
overdamped case without hydrodynamic interactions, the only important mode is the cage 
diffusion mode, i.e. the Fourier transform of the single-particle density fluctuations: 

N 

n(k) = ]T (e* k r * - <e lk - r %) (B.l) 
i=i 

In this case the lowest order mmca takes into account bilinear products of cage diffusion 
modes: n(k)n(— k)I 44 l. 

We start from the Green-Kubo expression eq.(34) for r/(0,a>) and eq.(35) for the stress- 
stress autocorrelation function p v (t;(f)). The first approximation of the mmca corresponds 
to the replacement of the full evolution operator e by its projection onto the subspace of 
the product variables n(k)n(— k) 



Here Q, is the N-particle Smoluchowski operator (cf. eqs.(35) and (37)) and P the normalised 
projector operator defined by 



where S eq (k;4>) = -^(n(k)n(- k)) e? is the equilibrium static structure factor and k runs 
over the reciprocal lattice. From eqs.(35), (B.2) and (B.3) we find for the stress-stress 
autocorrelation function 



atomic liquids! 19 !. The basic idea behind the mmca is that fluctuations (or 'excitations') 



(B.2) 




(B.3) 



Pr,(t;<i>) = 

k,k' 



(Sg,n(k)n(-k)) eg (n(k)n(-k) e m n(kOn(-kO) e9 (n(kOn(-kQSg ;; ) 
AN±S%{k-A)S 2 eq {k'-<P) 



(B.4) 
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The second approximation is to assume that the two modes appearing in the product 
variables propagate independently from each other. This means that the four-variable cor- 
relation function (n(k)n(— k)e nt n(k')n(— k')) eq in eq.(B.4) can be factorised into products 
of two- variable correlation functions (as already used in the normalisation of P (eq.(B.3)) 
giving 

(n(k)n(-k)e nt n(k')n(-k')) eq = 
= (n(k)e m n(-k')) eq (n(-k)e m n(k')) eq + (n(k)e m n(k')) eq (n(-k)e nt n(-k')) eq = 
= iV 2 F e 2 g (k;t)(5 kik , +(5 ki _ k (B.5) 

with F eq (k;t) = j^(n(k)e nt n(—k)) eq the equilibrium intermediate scattering function. As 
outlined in section 2 the main diffusion process at long times and high concentrations 
0.3 < 4> < 0.55 is the cage diffusion process, rather than free diffusion. Thus the long time 
decay of the equilibrium intermediate scattering function is determined by (&:;</>), the 
lowest eigenvalue, given by eqs.(13) and (14), corresponding to the eigenfunction n(k) of a 
kinetic operator defined elsewhere^ 11-13 !. This gives 

F eq (k; t) = S eq (k; 0) e —"( fc ^ (B.6) 

Performing the summation over k' and changing the summation over k to an integral over 
k in the limit of large volume V, we find from eqs.(B.4)-(B.6): 

167T3 J [ S eq {kAY 

where 



V„(k) = i < E3,n(k)n(-k) >„ (B.8) 

is the strength of the coupling between the microscopic stress tensor Y^. y (eq.(36)) and two 
microscopic densities. To evaluate V(k) we use that for an arbitrary function f(r N ) one 
has: 

< Fif(r N ) > eq = -k B T < Vif{r N ) > eq (B.9) 
32 



where r N = r±, - ■ ■ ,rjy. Eq.(B.9) follows from partial integration and using the explicit 
form of the equilibrium distribution function. Substituting eq.(36) for J2xy i n (B-8) and 
using (B.9) yields 

k T N r) 
V v (k) = — B — £ < ri , x _„(k)n(-k) > eg (B.10) 

iV i=i ari .y 

From (B.l) for n(k) and the expression below (B.3) for S eq (k; cf>) follows straightforwardly 



V v (k) = -k B Tk y -^-S eq (k; <j>) (B.ll) 



or equivalently, 



V v (k) = -k B T^S' eq (k; ( f ) ) (B.12) 

Substitution in (B.7) and performing angular integrations in k-space, leads to the final 
result for p v (t;4>), i.e., 

Then eq.(20) for r](4>,uj) follows immediately from eqs.(34) and (B.13). 
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Figure Captions 

1. Reduced cage-diffusion time t c (k; <p)/rp as a function of n = ka for volume fractions (j) 
= 0.30 (dotted line); 0,45 (dashed line); 0.50 (solid line) and 0.55 (dash-dotted line). For 
k ss k* ss 2tt the two times are of the same order of magnitude. 

2. Relative infinite frequency viscosity i] 00 (4>)/r]o as a function of the volume fraction (f>. □ 
Zhu et al (ref.22); x van der Werff et al (ref.5); • Cichocki and Felderhof (ref.23) whose 
points were obtained by a different analyses of van der Werff et al's than by the authors 
themselves (cf. Table II). The solid line corresponds to eq.(24). 

3. Relative Newtonian viscosity t}n(<I>)/t)o as a function of the volume fraction <f>. x van 
der Werff and de Kruif (ref.6); A van der Werff et al (ref.5) (cf. Table II); • Jones et al 
(ref.40); □ Papir and Krieger (ref.41). The solid line corresponds to eq.(25) and the dashed 
line to Voo{<t>)/rio = x{4>) (eq.(24)). 

4. Real (a) and imaginary part (b) of the reduced viscosity r} R ((j),uj) resp. rfj{(f),uj) as a 
function of 0JT\((j)). Experimental points from van der Werff et al (ref.5), © for <j> = 0.44, 
o for (j) = 0.46, □ for = 0.47, □ for <f> = 0.48, \/ for (f) = 0.51, * for 4> = 0.52, x for <f> = 
0.54 and A for <p = 0.57. Theory from eqs.(20), (25) and (30). Dashed line: <p = 0.55; solid 
line: <j) = 0.50; dotted line (j) = 0.45. The cloud of points in (b) near ujti(<P) = 1 should be 
discarded since they do not satisfy the Kramers-Kronig relation (ref.23). 

5. Relative real and imaginary parts of the visco-elastic viscosity, respectively: rf((/>, uj)/rjo 
(o) and rj"((f), uj)/rfo (x), as a function of ojt\{4>), for eight suspensions studied experimentally 
by van der Werff et al (ref.5) from 4> = 0.44 up to = 0.57 (cf. Table II). In order to make 
a fair and realistic comparison of the theory with experiment, keeping in mind the 4% 
uncertainty in the determination of <p and the extreme sensitivity of the denominator of 
r]* RI {(f),uj) - as already pointed out by van der Werff et al' 5 ! - we assign to the experimental 
data an effective volume fraction <p* , such that (r] t ^ eory '(</>*) — rf^ eory '{<p*)) = (^ P (0) — 
if^ v ((f))) , within the experimental uncertainty of (f>. Dotted line: phenomenological results 
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by Cichocki and Felderhof (ref.23) (only available for cp = 0.46, 0.54 and 0.57); solid line: 
theory from eqs.(20) and (29) using <(> = (/>* (cf. Table II). 

6. Ratio of T\{4>) and rp as a function of the volume fraction <\>. Experimental points from 
van der Werff et al (ref.5) (cf. Table II). Dashed line: theory from eq.(33); solid line: theory 
using eq.(41) instead of eq.(32) in eq.(30) in order to get the correct coefficient of the square 
root singularity at large frequencies (cf.section 6 and fig. 7). 

7. Coefficient of the square root singularity at large frequencies A(cp) as a function of the 
volume fraction <p. Experimental points from van der Werff et al (ref.5) (cf. Table II). 
Dashed line: mode-mode coupling theory (eq.(32)); solid line: exact result starting from 
the Green-Kubo relation (eq.(41)); dotted line: the theoretical result with Dq instead of 
D s (4>) = Dq/x{4>) (cf.section 6) 

8. (a) Inverse relative infinite frequency viscosity ^oA/ooO^) (• experimental points from 
van der Werff et al (ref.5); dashed line: theory from eq.(24)) and inverse relative Newtonian 
viscosity r/o/??iv(</>) (x experimental points from van der Werff et al (refs.5 and 6); solid 
line: theory from eq.(25)) as a function of the volume fraction <p. Dotted line: Beenakker's 
expression (44c) (ref.30) (cf. Section 8, sub 3(a)). 

(b) Relative short time self-diffusion coefficient D s (^>)/Dq as a function of the volume 
fraction cf>. □ Zhu et al (ref.22); x Van Megen et al (ref.45); • Pusey and Van Megen 
(ref.46). The solid line corresponds to eq.(45a) and the dashed line to the Beenakker and 
Mazur expression (45b) (ref.31). 

(c) Inverse relative infinite frequency viscosity t/oA?oo(<A) (• Zhu et al (ref.22); □ Van 
der Werff et al (ref.5)) and relative short time self-diffusion coefficient D s (4>)/Dc) (o Zhu et 
al (ref.22); □ Van Megen et al (ref.45)) as a function of the volume fraction <p. Solid line: 
theory from eq.(47b); dotted line: Beenakker (ref.30); dashed line: Beenakker and Mazur 
(ref.31). 
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Table I 

Characteristic values of the model systems used! 5,6 ] 



System 


a(nm)(DLS) 


r^lO^s-W 2 ) 


Tp (ms) 


SP 23 


28 ± 2 


8.68 


0.0903 


SSF 1 


46 ± 2 


5.29 


0.400 


SJ 18 


76 ± 2 


3.20 


1.81 



Table II 
Parameters discussed in text. 





System 


ti{4>)/t p 


r)oo{<t>)/Vo 


vn(4>)/vo 


A(</>) 


<t>* 


0.44 


SSF 1 


0.402 


4.99 


12.2 


7.69 


0.431 


0.46 


SP 23 


0.421 


5.13 


13.1 


8.33 


0.438 


0.47 


SJ 18 


0.776 


6.78 


17.8 


8.45 


0.458 


0.48 


SSF 1 


0.372 


6.36 


17.3 


12.1 


0.458 


0.51 


SJ 18 


0.665 


7.45 


28.8 


17.7 


0.498 


0.52 


SSF 1 


0.834 


7.47 


32.7 


18.6 


0.508 


0.54 


SSF 1 


0.912 


9.9 


50.7 


28.8 


0.535 


0.57 


SSF 1 


3.70 


11.5 


139 


44.7 


0.593 


0.58 


SP 23 


3.99 


10.0 


187 


60.2 
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Figure 5a-d 
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Figure 5e-h 
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Figure 8a 
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Figure 8b 
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